Universal metallic and insulating properties 
of one dimensional Anderson Localization : 
a numerical Landauer study 
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We present results on the Anderson localization in a quasi one-dimensional metallic wire in the 
presence of magnetic impurities. We focus within the same numerical analysis on both the uni- 
versal localized and metallic regimes, and we study the evolution of these universal properties as 
the strength of the magnetic disorder is varied. For this purpose, we use a numerical Landauer 
approach, and derive the scattering matrix of the wire from electron's Green's function obtained 
from a recursive algorithm. 



Interplay between disorder and quantum interferences 
leads to one of the most remarkable phenomenon in con- 
densed matter : the Anderson localization of waves. The 
possibility to probe directly the properties of this local- 
ization with cold atoms [H [3] have greatly renewed the 
interest on this fascinating physics. In this paper, moti- 
vated by the transport properties of metallic spin glass 
wires 0, 0, Q, we study the electronic localization in 
the presence of both a usual scalar random potential and 
frozen random magnetic moments. 

One dimensional disordered electronic systems are al- 
ways localized. Following the scaling theory [6| this im- 
plies that by increasing the length L x of the wire for a 
fixed amplitude of disorder, its typical conductance ul- 
timately reaches vanishingly small values. The localiza- 
tion length £ separates metallic regime for small length 
L x <C £ from the asymptotic insulating regime. In the 
present paper, we focus on several universal properties 
of both metallic and insulating regime of these wires in 
the simultaneous presence of two kind of disorder, orig- 
inating for example from frozen magnetic impurities in 
a metal. The first type corresponds to scalar potentials 
induced by the impurities, for which the system has time 
reversal symmetry (TRS) and spin rotation degeneracy. 
In this class the Hamiltonian belongs to the so-called 
Gaussian Orthogonal Ensemble (GOE) of the Random 
Matrix Theory classification (RMT). If impurities do 
have a spin, the TRS is broken as well as spin rotation 
invariance. The Hamiltonian is then a unitary matrix, 
which corresponds in RMT to the Gaussian Unitary En- 
semble (GUE) with the breaking of Kramers degeneracy 
Q. However, for the experimentally relevant case of a 
magnetic potential weaker than the scalar potential, the 
system is neither described by the GUE class, nor by the 
GOE class, but extrapolates in between. 

In this paper, we study numerically the scaling 
of transport properties of wires with various relative 
strength of these two disorders. The phase coherence 
length which phenomenologically accounts for inelas- 
tic scattering^ is assumed larger than the wire's length 



L x . We describe the disordered wire using a tight-binding 
Anderson lattice model with two kinds of disorder poten- 
tials : 
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The first scalar disorder potential V — {v.i}i is diagonal 
in electron-spin space. Vi is a random number uniformly 
distributed in the interval [— W/2, W/2]. s,s' label the 
SU(2) spin of electrons and the Si correspond to the 
frozen classical spin of impurities, with random orien- 
tations uncorrelated from an impurity site to another. 
Varying the amplitude of magnetic disorder J allows to 
extrapolate from GOE to GUE. For a given realization 
of disorder, the Landauer conductance of this model on 
a 2D lattice of size (in units of lattice spacing) L x x L y 
is evaluated numerically by a recursive Green's function 
technique. We stay near the band center, avoiding the 
presence of fluctuating states studied in [lOJ . Universal 
properties are identified by varying the transverse length 
L y from 10 to 80, with the aspect ratio L x /L y taken from 
1 to 6000. Typical number of disorder averages is 5000, 
but for L y — 10 we sampled the conductance distribution 
for 50000 different configurations of disorder. 

Localization length. The localization length £ is ex- 
tracted from the scaling of the typical conductance in 
the localized regime according to [111, [l2| : 



exp(logg) = exp {-2L X /C) . 



(2) 



where g = G/(e 2 /h) is the dimensionless conductance 
and (•) represents the average over scalar disorder V. Nu- 
merical plot of the logarithm of this equation is given in 
the inset of figured] The linear behavior of (log (?) is high- 
lighted in the insulating regime and fitted to provide the 
localization length for each value of the transverse length 
L y . Note that we checked that the very slow convergence 

of the Lyapunov exponent j(L x ) 



g{L x ) J 
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towards £ 1 {L y ) provided comparable results. The evo- 
lution of £ with the transverse length is expected to follow 



£ = (PL y + 2 - 0)l e 



(3) 



with l e the mean free path and (3 = 1 corresponds to 
the orthogonal universality class GOE while (3 = 2 for 
GUE. Note that this change in (3 is accompanied by 
an artificial doubling of the number of transverse modes 
N y = L v — > 2N y due to the breaking of Kramers degen- 
eracy [111 ]. Comparison of numerical localization lengths 
for different J with is shown in fig[TJ Excellent agree- 
ment is found for J = (GOE class, (3 — \). In the case 
J ^ Owe observe a crossover between GOE and GUE for 
intermediate values of magnetic disorder, while a good 
agreement with the GUE class is reached for J > 0.2. 
From these results, we already notice that the localiza- 
tion regime is reached for much longer wires in the GUE 
case than for GOE. As shown below, this allows for an 
easier numerical investigation of the universal metallic 
regime in the GUE case : magnetic impurities help in 
finding the universal conductance fluctuations ! 




FIG. 1: Evolution of localization length as a function of trans- 
verse length. l e is the mean free path of the diffusive sample. 
Different behavior of the localization length if J = or J 7^ 0. 
Inset : Scaling of the typical conductance (log g) = — 24*. 



Insulating regime. In the insulating regime L x > £, 
we expect a Log-normal conductance statistical distribu- 
tion [llj |. However in the region g = 1 and for (g) < 1, we 
find a non-analytical behavior of P(g) in agreement with 



13l Il4 lla . 1 16l ] as shown for instance in fig. [51 In this 



figure we plot the distribution P{g) for similar values of 
(g) for GOE (J = 0) and GUE (J = 0.2). The shapes of 
these distributions are highly similar if (g) <C 1, showing 
that both distributions tend to become Log-normal with 
the same cumulants. In the intermediate regime, shapes 
are symmetry dependent. Moreover the non-analyticity 
appears for different values of conductance (close to 1) 
and the rate of the exponential decay [l3| in the metallic 
regime seems to differ from one ensemble to the other 



<g) = o.6S 




FIG. 2: TOP: Comparison of Probability density functions 
(PDF) of conductance for J — (plain curves) and J = 0.2 
(dashed curves). Plots are performed for different values of 
average conductance, (a): {g)(J = 0) = 0.84 and (g) (J = 
0.2) = 0.79. (b): (g)(J = 0) = 0.67 and (g)(J = 0.2) = 
0.62. (c): {g)(J = 0) = 0.45 and (g)(J = 0.2) = 0.42. (d): 
(g)(J = 0) = 0.21 and (g)(J = 0.2) = 0.18. BOTTOM: PDF 
of conductance for (g) < 1(J = 0.2) and (g) > 1(J = 0) and 



Gaussian interpolations. L y 
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(see for instance curves (a) or (b)). Finally, the bottom 
curve of figure represents the distribution of conduc- 
tance for just above and below the threshold (g) = 1. 
Plain lines represent gaussian interpolations with a mean 
and a variance given by the first and the second cumulant 
of each numerical conductance distribution. For (g) > 1, 
the gaussian interpolation approximates very well the full 
distribution. On the other hand, as soon as (g) < 1, the 
gaussian law only approximates the tail g > 1 of the 
distribution of conductance . This behavior is in agree- 
ment with the sudden appearance of the non-analyticity 
for distributions with average conductance inferior to 1 
[l6[ . This conductance distribution converges to the Log- 
normal only deep in the insulating regime, the conver- 
gence being very slow (much slower than in the metallic 
regime) . This qualitative result is confirmed by the study 
of moments : in the insulating regime the second cumu- 
lant is expected to follow [17j: 



((log ff -(log 3 )) 2 ) H(log 5 ) 5 



-2(log.g), (4) 



FIG. 3: Plot of the variance of log g as a function of the mean 
for the orthogonal and unitary case. In inset is shown the 
third cumulant of log;; as a function of the mean. 



Our numerical results are in agreement with this scaling 
(figure [3]) with however very slow convergence towards 
this law : corrections are measurable even if the system is 
deeply in the localized state. More precisely, we find (see 
fig. [3]) that for the deep insulating regime ((logg) 2 ) c = 
— 1.88(log<?) slope —1.88, with a slight discrepancy with 
. Finally in the inset of figure [3J we show the third 
cumulant of logg as a function of the first one. The 
linear behavior is in agreement with the single parameter 
scaling. We find that contrary to the second cumulant 
the coefficient of proportionality between the skewness 
and the average depends on the symmetry of disorder, 
which is not expected. 

Metallic regime. We now focus on the universal 
metallic regime described by weak localization. By defi- 
nition weak localization corresponds to metallic diffusion, 
expected for lengths of wire l e <C L x <C £. For this regime 
to be reached, we thus need to increase the number of 
transverse modes L y and thus £ for all other parame- 
ters fixed (see (|3])). Moreover, for a fixed geometry, this 
regime will be easier to reach in the GUE class than in 
the GOE. Fig. [2] shows that the conductance distribu- 
tion is in good approximation Gaussian with a variance 
described [9( by 




(5) 



where x = L x /L m and the scaling function F{x) depends 
only on dimension [{j fl8|. In the bottom plot of figure |H 
these conductance fluctuations are plotted as a function 
of longitudinal length L x for different values of J. A 
single parameter fit by (|5|) provides the determination of 
the magnetic dephasing length L m as a function of the 
magnetic disorder J. The determination of L m allows 
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FIG. 4: TOP : second cumulant of g as a function of first 
cumulant of g showing the universal behavior in the metallic 
regime. BOTTOM : variance of g as a function of longitudinal 
size. UCF are shown. Different curves correspond to different 
values of magnetic disorder J. In inset variance of g is plotted 
as a function of L x /L m . Transverse length L y — 40. 



for very fine comparison with weak localization theory 
in this regime, allowing for example the study of con- 
ductance correlation between different disorder configu- 
rations (see 18,1||). The inset of Fig. 2]shows the scaling 
form of these fluctuations (as a function of L x / L m (J)) 
in excellent agreement with the theory Moreover, 



for long wires (and large values of J) conductance fluc- 
tuations are no longer L x dependent and equal to 1/15. 
This is the so-called Universal Conductance Fluctuations 
(UCF) regime which is precisely identified numerically in 
the present work. The top plot of figure 2] confirms an- 
alytical results from [l5| both qualitatively in the shape 
of the curves and quantitatively in the values of fluctua- 
tions in both universality classes. In our study, values of 
UCF are reached with a maximal error of 1% for GOE 
and 3% for GUE with respect to the analytical value 
of the UCF in the regime independent of (g) (i.e with 
much higher precision than e.g [20( and 14 1) . The other 



information provided by this curve is the condition for 
having a universal behavior, i.e independence on the ge- 
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FIG. 5: Plot of (g 3 ) c as a function of (g) in the metallic 
regime. In insets are shown convergence curves for L y = 10 
andJ = or J = 0.2. 

ometry or on the disorder of the sample. As the localized 
regime is harder to reach for GUE, the universal metallic 
regime is easier to obtain, the plateau of UCF is widen. 
One needs then larger transverse lengths (recall that £ 
increases with L y ) to get the UCF regime for GOE, this 
is exactly what we see on figure 2] where the UCF plateau 
is reached only for L y = 80 for GOE whereas it is reach 
for L y = 80 and L y = 40 for GUE. 

Finally we consider the third cumulant of the distribu- 
tion of conductance. According to the analytical study 
of [151 ] , this cumulant decays to zero in a universal way as 
(g) increases. Here we find a dependance of this decrease 
on the symmetry class : for GOE (<? 3 ) c goes to zero in 
a monotonous way whereas it decreases, changes its sign 
and then goes to zero in GUE case. For (g) > 4 numer- 
ical errors are dominant, then this part of the curve is 
irrelevant. Note that this fast vanishing of the third cu- 
mulant confirms the faster convergence of the whole dis- 
tribution towards the gaussian, compared to what hap- 
pens in the insulating regime. Based on our numerical 
results, we cannot confirm nor refute the expected law 
i£) c oc l/(g} n , with n = 2 in GOE and n = 3 in GUE 

To conclude we have conducted extensive numerical 
studies of electronic transport in the presence of ran- 
dom frozen magnetic moments. Comparing and extend- 
ing previous analytical and numerical studies, we have 
identified the insulating and metallic regimes described 
by the universality classes GOE and GUE. We have paid 
special attention to the dependance on this symmetry 
of cumulants of the distribution of conductance in both 
metallic and insulating universal regimes. In particu- 
lar, we have identified with high accuracy the domain of 
universal conductance fluctuations, and determined its 
extension in the present model. 
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Note : during the very last step of completion of this 
paper, we became aware of a preprint by Z. Qiao et al. 
[231 ] which performed a similar numerical Landauer study 
of ID transport for various universality classes, and fo- 
cused mostly on the metallic regime. While both studies 
agree on the finding of UCF (although we have higher 
accuracy for (3 — 1), we did not find signs of a second 
universal plateau for ((Sg) 2 ) in our study. 
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